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This research used images from Optical Coherence Tomography (OCT) 
examination as well as fundus images to localize the optical disc and macular 
layer of retina. The researchers utilized the OCT and fundus image to 
interpret the distance between macular center and optic disc in the image. 
The distance will express the area of macula that can be employed for further 
research. This distance could recognize the thickness of macula parameters 
diameter that will be used in localizing process of optic disc and macula. The 


Keyword: parameters are the circle radius, the size of window’s filter, the constant 
value and the size of optic disc element structure as well as the size of 
Fundus . i 
macula. The results of this study are expected to improve the accuracy of 
Macula macula detection that experience the edema. 
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1. INTRODUCTION 

One of the indicators for diabetic retinopathy in advance stage is edema macula. It is a leakage of 
fluids in retina. Diabetes can cause small blood vessel in this area to release fluid or fat. Edema macula may 
affect the patient’s vision center and the worst-case edema macula may cause blindness [1]. Macula is a part 
of retina that allows eye to see soft details with a spherical morphology that located in the back of retina at 
3.5 mm temporal distance and 0.5 mm smaller to the discus [2]. In the macula, there are vessels that play 
important role as photoreceptor. The study in sharpening the retinal fundus images was performed by the 
researchers using adaptive histogram and morphological operation of tophat as well as bothat. Blood vessel 
segmentation is necessary to perform in localizing other objects of retinal fundus images. These images are 
optic disc and macula [3]. In addition, retinal vascular extraction also performed using curvelet 
transformation and morphological operation that combine two operations. Those operations are tophat dan 
bothat [4]. The localization of the optic disc and macula area must be seperated from the blood vessel. 
Therefore, the sharpness of human vison is highly dependent on this section [5]. The macula lies in the mid- 
posterior of retina. 

One of the initial phase in performing macula detection that experienced edema is by localizing the 
area of optic disc and macula area of retinal image. Optic disc localization in this study was carried on 
through tagging the connected component based on color intensity that have similarities with exudates by 
checking pixel (P) per pixel (from left to right side and from top to bottom direction) to identify the 
connected pixels area. The connected pixel area is an area from the bordered pixel that have the same 
intensity or intensity value in a set of V (in binary image V = {1}). The connected component for input 
image follows the connectivity rules (4, 8 or m-connectivity) [6]. This study proposes a new technique to 
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analyze the swelling of the macula. This study consists of two stages. The first stage is vessel extraction and 
the second one is vessel reconstruction. 

The image input is a fundus green channel image that is converted into grayscale image. In this 
study, Speeded Up Robust Features (SURF) algorithm was used to describe the process of macula detection 
on retinal image by discovering the compatibility between the features points of each image. This study also 
used Random Sample and Consensus (RANSAC) algorithm to match the points that correspondent and 
Gaussian function algorithm for smoothing the macular image splicing. The RANSAC algorithm that was 
implemented in this study used to select random subset from original image to achieve the best results in 
localizing the swollen macula area [7]. This study determined the presence of macula edema by detecting the 
presence of fovea. Fovea lies in the middle of macula. Fovea and optic disc detection were performed using 
motion pattern recognition. The feature selection was carried on using Principal Component Analysis (PCA) 
by creating new feature derived from extraction of linear features combinations. Therefore, they can be used 
to reduce the selected feature dimension which is edema macula. Gaussian derivative also used in this study 
to capture the characteristic of edema macula patters that will be detected [8]. 

The retinal fundus optic disc locatization was performed in several steps. The first step was 
performing green chanel image fundus extraction. The next step was carry out fitur optic nerve head (ONH). 
ONH was performed based on the contour from optic disc cup area. In performing the extraction, the 
researchers converting the image to grayscale and executing thresholding with dilation and erotion operation. 
The next step is performing edge detection using dilation operation to create matrix feature. The result of 
matrix feature will be used as an input for training and testing process. The training and testing process in this 
study using SVM method [9]. The research was performed by establishing an automatic diagnostic system to 
detect macule using fluorescein image. In the preliminary stages, noise reduction was done in preprocessing 
phase to eliminate the object details beside the macula and to maintain its lightning distribution with 
Gaussian filter. There were several processes in this study, for instance blood vessel segmentation with 
Gaussian filter, optic disc detection, and red lesion detection to segment micro aneurysm feature. 

The segmented feature will describe the shape, color and contrast of the detected candidates [10]. 
The optic disc detection in this study was done by detecting the maximum area with the highest pixel value. 
Then the retinal fundus image along with the optic disc that already being detected will perform local image 
features process. Pixel weighting was carried out to gain neighboring influence in performing segmentation 
by clustering each pixel inside the cluster center. The macula classification phase proposed in this study is 
based on some parameters such as the size of macula window’s area and spatial relationship between pixels 
in the retinal fundus [11]. 

The medical examination of the edematous macules is done by direct observation by the physician 
through Optical Coherence Tomography (OCT) examination. This examination will provide cross sectional 
of retinal issue images. In this examination, the thickness of the retina and macula will be measured, which 
indicate the occurrence of swelling and fluid clots in the area. Determination of optic disc localization in this 
study carried on using fundus image. Disc localization and macula are principal factors that determined the 
distance interpretations between the centers of macula with optic disc in pixel representation to know the 
diameter of macular thickening. This study used some parameters. For instance, the circle radius, the size of 
window filter, and constant value as well as the size of optic disc structuring element and macula. The result 
of this study is expected to improve the accuracy of macular detection that experienced edema. 


2. RESEARCH METHOD 
2.1. Data Sets 

This study used imaging images with Optical Coherence Tomography (OCT) and Fluorescein 
Angiography fundus image as can be seen in Figure 1. 


Figure 1. OCT and Fundus Image 
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In the fundus retinal image in Figure 1, the image shows quite low contrast of retina’s macula. 
Macula objects are sometime hindered by another object that indicate retinal abnormalities such as the 
exudate, dot and blood hemorrhage, micro aneurysm as well as Intra-Retinal Microvascular Abnormalities 
(IRMA’s). These abnormalities can cause difficulties and in localization of macular area. 


2.2. Estimate Optic Disc Position 
In Fluorescein Angiography (FA) fundus image, the initial phases of the process were detecting the 
center estimation and the diameter of optic disc. The estimation of macula area was based on the distance and 
its position with optic disc. This distance is relatively constant in each retina. The size of the disc on each 
retina image of the patient is varies, but generally the size of optic disc has diameter between 1.8 + 2mm 
inside the image fundus with standard size of 768 x 576 pixel. There are some factors that make optic disc 
edge detection accurate. Some parts of the optic disc will be obscured by many retinal blood vessels. In 
general, the process in localizing the existence of the optic disc in this study was performed by the following 
algorithm: 
1. Green channel FA fundus image extraction. The green channel was selected as the primary color since 
green color has the high intensity compare to the others two colors; red and blue in an image with RGB) 
color as can be seen in Figure 2. Below is the mathematical formula to find the green channel selection : 


G () 


£= (R+G+B) 


Figure 2. OCT and Green Channel Fundus Fluorescein Angiography (FA) Image 


2.3. Morfology 

Formula (1) is performed to find a green channel in an image, while R, G and B are stand for Red, 
Green and Blue sequentially. In green channel, the detail of image can be seen clearly and comprehensively. 
Whereas the used of red channel only going to show the image edge while in blue channel merely provide 
just part of the image’s detail and contain a lot of noise. Based on image in Figure 2, the contrast of optic disc 
can be observed better in green channel compared to red and blue channel in RGB image. Therefore, this 
study only used green channel image to compute the optimal threshold. Figure 3 shows the green channel 
fundus Fluorescein Angiography (FA) image and its histogram. It can be seen that the pixel that 
corresponding with optic disc have higher intensity in histogram. 
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Optic Disc 
Arsa Histogram 


Figure 3. OCT and Green Channel Fundus Fluorescein Angiography (FA) Image 


2. The localization of optic disc position. Optic disc localization can be performed by considering the 
geometry shape of optic disc on several connected components. The difficulties in estimating the optic 
disc position is the presence of other objects outside the optic disc image that has connected component 
as well as can be seen in Figure 4. 


Optic Dise exudates 


Figure 4. Connected Component Detection on Optic Disc 


3. Each component that connected to the image will be labelled with the number of pixels in the 
component. Spatial coordinate of each connected component is calculated. Component is assumed to 
have maximum number of pixels with maximum diameter of optic disc which is 2 mm. Consequently, if 
one component of the image spatial coordinate is within 50 distance pixels from spatial coordinate of the 
mean largest component, the object will be combined and the new mean of image spatial coordinate will 
be computed as can be seen in Figure 5. 


Figure 5. Position Estimation of Optic Disc 


4. Performing connected component tagging with object geometry approach using some parameters below : 
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1) The value area is computed based on the maximum and mean intensity in the optic disc area. 

2) Restricting the area that are detected as an optic disc. This area is the number of image dimension 
with angle and width parameters in retinal image that is written in an array form using the formula: 
BoundingBox=[corner width] (2) 

3) Retinal optic disc position estimation based on the ratio where the shape of the optic disc is circle. In 
a circle object the properties that can be discovered are [12]: Radius of the circle r , circle diameter 
d , wide of circle L , circumference of circle and P parameter. The formula to localize optic disc 
position based on the ratio of round 

4) The extensive and parameter value as the property of this circle can be performed in optic disc 
regions which will be localized as the basic form of the roundabout size using this formula: 

Ratio = 2— 3) 
7a 
5) If the result value of detected object ratio yields 1 and has the intensity equal or almost the same as 


optic disc intensity, this detected object will be marked (localization). The objects that have different 
characteristic with the optic disc will not be marked. In general, the optic disc illustration can be 
seen in Figure 6. 


Figure 6. Optic Disc Position Localization Illustration 


The ratio value of | is based on the ratio length of major axis and minor axis of the object that was 


detected in fundus image as can be seen in Figure 7. 


Figure 7. Ma Detected Optic Disc 
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2.4. Localization OCT Macula Layers 
Localization of macular image layer area was performed on preprocessing phase. This localization 
was carried on in calculating the macular area outside the normal size of retina macular image. The macula 
area has the thickness of 400 um, diminishing on fovea with the size of 150 um, and becoming thinner on ora 
serrata with the thickness of 80 um [13]. The macula has very low contrast inside the retinal image and 
sometime the area of macula is blocked by the exudates or the bleeding in the retina’s region [14]. The 
algorithms of preprocessing phase on OCT macular images area can be explained below: 
1. Read the OCT Image 
2. Performing the blurring process to reduce the noise in OCT image. The intensity value of the image 
changes a lot as in Figure 8. The noise in OCT images have the form of dots or pixels with different 
intensity. This study used the first derivative of Gaussian filter to reduce the noise on OCT image with 
the formula below: 


G(x, y)= ze 20° (4) 


Where Ø denote the derivative standard from pixel distribution of OCT image OCT Į (x, y). Assumed a 


distribution have an average value equal to 0 ( O =1). 


Figure 8. Noise Reduction Process and OCT Macula Image 
(a) Original Image (b) Image from Blurring Process 


3. Performing edge detection on image gradient at horizontal direction from the filtering image using the 
algorithms below: 
a) Read the image resulted from blurring phase. 
b) Read size of the image with function [N, M] =size (G). [N M] is a pixel line and column of an image 
obtained from blurring phase G(x,y). The first derivative of image function G(x,y) towards xand y 
axis that are computed using formula [15] : 


8G ex, y) and oes. y (5) 
ox oy 


Whereas, the vector gradient and the amplitude can be calculated with [9] : 


2G x,y) 
VG(x,y) =| & o 
sae (x, y) 
oy 
|VG|(x, = (Zh yy + (Zh y)? 7) 
Ox oy 


4. Referring to equation 7, the first partial derivative of images from blurring phase were computed based 
on the two-closed neighboring pixel. Therefore, the gradient can be obtained: 
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VG(x, y), = < (x, y)=G(x+1, y)- G(x, y) (8) 


VG(x, y), is an image gradient on x horizontal direction that produced edge of object in the shape of 


horizontal line, as in Figure 9. 


Figure 9. Image Gradient on Horizontal Direction 


5. Labeled the connected component that link horizontal image gradient with adjacent matrix. x pixel is 
said to be adjacent with y pixel if there is a connected connection between them. Two images subset Sı 
and S2 are adjacent if some pixels on S; adjacent with some Sı pixels as in Figure 10. 


Figure 10. Subset Sı and S2 


Path from pixel x with coordinate (p,q) for y pixel with coordinate (7, 5) is a pixel order that differ from 


coordinate (p,,qy),(Py+41)s-(Pn In) 
where: 


n = Path 
(P0240) =(P1 Ds (Pio In) = (1,8) (9) 


(p,,y,) is adjacent where : 
Piao fia (10) 


6. Compute the distance between x( p,q) and y(r,s) using city block distance with 4 connectivity as 
can be seen in Figure 11 using the formula : 


D, xy) p-rl|+lq—-s| (11) 
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y 


Figure 11. City Block Distance with 4 Connectivity 


3. RESULTS AND ANALYSIS 
3.1. Experimental Image 

The localization of the optical disc position was performed in geometric optic disc form by 
considering several connected components. The difficulty in performing the optic disc position estimation 
was the occurrence of another object outside the optical disc that has a connected component as well. Figure 
12 shows the proposed algorithm in this study succeed to localize the presence of optic disc on fundus image. 


Non-Optic Disc 


Optic Disc 


Figure 12. Estimation Image OD Position 


In previous study, the process of retinal fundus optic disc images localization was performed using 
combination of morphological operation and optic disc center point determination. This method was carried 
out by detecting the geometric shape in the parametric form, where the shape of the optic disc is circular. The 
exploitation of Hough Transformation [16] was carried out to search the most suitable geometric form that 
correspond to the set of points in the image. The transformation is done to find the lines and the curves that 
passing by many points inside the image. These lines and curves are the one that near and most correspond 
with the data in the image matrix. The selection of the circular radius greatly affect the accuracy of the area 
localization of optic disc. If the circular parameter value is too small or too big, it will affect the distance of 
macula’s center from the central of optic disc, resulting macula ROI will be unsuitable. The study in this 
paper used roundness ratio. This roundness ratio is the ratio of the radius at the smallest angle to the radius 
average of fundus image as object. The ratio of this study was proposed by the researchers to perform 
comparison between the shorthest axis with the longest axis of the optic disc. Table 1 shows examples of FA 
fundus Image OD estimation position. 

As can be seen on Table 1, the optic disc with the roundness object ratio approach succesfully 
localizing its position correctly. This ratio is the average geometric diameter of optic disc with its longest 
diameter, where the object with a rounded geometric shape connected component in fundus image is an optic 
disc. Therefore, those area are localized as optic disc. 
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Table 1. Optic Disc Localization Result 


Image FA Fundus Image Optic Disc Localization Result 
7 | 
7 [ 
i e 
7 EJ 


3.2. OCT Macula Image Localization 
The macula’s layers play important role in human’s vision. These multiple layers are located in 
three different areas (Fovea, Avascular Zone and Foveolar) as can be seen in Figure 13. 


Figure 13. (a) Fundus Scanned Image FA, (b). Area and Layers of OCT Macula, 
(c). Macula Area Position 


Edema macula is associated with vision blindness in non proliferative diabetic retinopathy. The 
symptoms are characterized by (i) the retinal thickening at or within 500um distance from the center of 
macula, (ii) the occurrence of exudates or within 500um distance from the center of macula followed by the 
thickening of the retinal surround the macula’s center and (iii) retinal thickening that greater than one disc 
diameter (DD) distance from the center of macula. The changes inside the macula’s layer area will influence 
the peripheral’s vision. Therefore, it is very important to localize the area of macula in detecting macula that 
has edema [17]. Figure 14 shows the example of OCT macula layers localization’s results on Inner Limiting 
Membrane (ILM) area and RPE OCT macula image. 
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Code Note 
ILM Inner Limiting Membrane 
NFL/ GCL Nerve Fiber Layer 
IPL/ INL Inner Plexiform Layer 
INL/ OPL Inner Nuclear Layer 
OPL/ ONL Outer Plexiform Layer 
OPL/ ONL Outer Nuclear Layer 
IS/ OS Inner Segment 
RPE Retinal Pigmented Epithellium 


Figure 14. Macula Area Localization and OCT Macula Layer [17] 


Many approaches had been done by some researchers in OCT macula images localization. This 
localization was carried out through selecting the pixel that will be used as seed point and by searching the 
similar neighboring pixels. The pixels area then expanded by adding similar pixel to the neighboring pixels. 
Unfortunatelly, researchers facing difficulties of performing segmentation on OCT makula area. The 
difficulties are on color similarities and on texture of each macula layer. Macula layers consist of 8 layers 
OCT. In this study, the researchers proposed the utilizing of adjacency matrix to present the pixels of macular 
layer. These layers are closed to the other macula’s layer. The adjacency was performed by presenting pixel 1 
which indicating if node i to node j has an edge and 0 otherwise (when the node has no edge). 


4. CONCLUSION 

The proposed algorithm succeeds to localize the existence of optic disc. This optic disc is based on 
the connected component that distinguished optic disc with the other object of the image fundus such as the 
exudates. The result of macula’s area localization on OCT image can compute the macula area outside the 
normal size of macula retina image. The approach with geometric optic disc ratio succeeds to localized the 
optic disc. The exploitation of adjacendy matrix able to distinguish every layer of OCT macula image despite 
their similarities and color neighboring. 

For the further research, the thickness of 7 layers of macula OCT image will be detected. This 
thickness will be used as reference to detect the sign of the presence of edema macula on OCT image. 
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